function [ VecContrib ] = ERdC_VecContrib_croise( donnees,elementCourant )

VecContrib = zeros(1,4) ;
L = donnees.poutre.L ;
F = donnees.F ;

a = donnees.forceLineiqueA ;
b = donnees.forceLineiqueB ;
c = donnees.forceLineiqueC ;

xi = elementCourant(1) ;
xj = elementCourant(2) ;

%VecContrib(1) = F - c*L - (b*L^2)/2 - (a*L^3)/3 + (c*xi)/2 + (3*b*xi^2)/10 + (a*xi^3)/15 + (c*xj)/2 + (2*b*xi*xj)/5 + (1/10)*a*xi^2*xj + (3*b*xj^2)/10 + (1/10)*a*xi*xj^2 + (a*xj^3)/15 ;
%VecContrib(2) = F*L - (c*L^2)/2 - (b*L^3)/3 - (a*L^4)/4 - F*xi + c*L*xi + (1/2)*b*L^2*xi + (1/3)*a*L^3*xi - (5*c*xi^2)/12 - (7*b*xi^3)/30 - (a*xi^4)/20 - (c*xi*xj)/6 - (2/15)*b*xi^2*xj - (1/30)*a*xi^3*xj + (c*xj^2)/12 - (1/30)*b*xi*xj^2 - (1/60)*a*xi^2*xj^2 + (b*xj^3)/15 + (a*xj^4)/60 ;
%VecContrib(3) = -F + c*L + (b*L^2)/2 + (a*L^3)/3 - (c*xi)/2 - (3*b*xi^2)/10 - (a*xi^3)/15 - (c*xj)/2 - (2*b*xi*xj)/5 - (1/10)*a*xi^2*xj - (3*b*xj^2)/10 - (1/10)*a*xi*xj^2 - (a*xj^3)/15 ;
%VecContrib(4) = (-F)*L + (c*L^2)/2 + (b*L^3)/3 + (a*L^4)/4 - (c*xi^2)/12 - (b*xi^3)/15 - (a*xi^4)/60 + F*xj - c*L*xj - (1/2)*b*L^2*xj - (1/3)*a*L^3*xj + (c*xi*xj)/6 + (1/30)*b*xi^2*xj + (5*c*xj^2)/12 + (2/15)*b*xi*xj^2 + (1/60)*a*xi^2*xj^2 + (7*b*xj^3)/30 + (1/30)*a*xi*xj^3 + (a*xj^4)/20 ;

% VecContrib(1) = -F - c*L - (b*L^2)/2 - (a*L^3)/3 + (c*xi)/2 + (3*b*xi^2)/10 + (a*xi^3)/15 + (c*xj)/2 + (2*b*xi*xj)/5 + (1/10)*a*xi^2*xj + (3*b*xj^2)/10 + (1/10)*a*xi*xj^2 + (a*xj^3)/15 ;
% 
% VecContrib(2) = (-F)*L - (c*L^2)/2 - (b*L^3)/3 - (a*L^4)/4 + F*xi + c*L*xi + (1/2)*b*L^2*xi + (1/3)*a*L^3*xi - (5*c*xi^2)/12 - (7*b*xi^3)/30 - (a*xi^4)/20 - (c*xi*xj)/6 - (2/15)*b*xi^2*xj - (1/30)*a*xi^3*xj + (c*xj^2)/12 - (1/30)*b*xi*xj^2 - (1/60)*a*xi^2*xj^2 + (b*xj^3)/15 + (a*xj^4)/60 ;
% 
% VecContrib(3) = F + c*L + (b*L^2)/2 + (a*L^3)/3 - (c*xi)/2 - (3*b*xi^2)/10 - (a*xi^3)/15 - (c*xj)/2 - (2*b*xi*xj)/5 - (1/10)*a*xi^2*xj - (3*b*xj^2)/10 - (1/10)*a*xi*xj^2 - (a*xj^3)/15 ;
% 
% VecContrib(4) = F*L + (c*L^2)/2 + (b*L^3)/3 + (a*L^4)/4 - (c*xi^2)/12 - (b*xi^3)/15 - (a*xi^4)/60 - F*xj - c*L*xj - (1/2)*b*L^2*xj - (1/3)*a*L^3*xj + (c*xi*xj)/6 + (1/30)*b*xi^2*xj + (5*c*xj^2)/12 + (2/15)*b*xi*xj^2 + (1/60)*a*xi^2*xj^2 + (7*b*xj^3)/30 + (1/30)*a*xi*xj^3 + (a*xj^4)/20 ;

VecContrib(1) = -F - c*L - (b*L^2)/2 - (a*L^3)/3 + (c*xi)/2 + (3*b*xi^2)/20 + (a*xi^3)/15 + (c*xj)/2 + (b*xi*xj)/5 + (1/10)*a*xi^2*xj + (3*b*xj^2)/20 + (1/10)*a*xi*xj^2 + (a*xj^3)/15 ;

VecContrib(2) = (-F)*L - (c*L^2)/2 - (b*L^3)/3 - (a*L^4)/4 + F*xi + c*L*xi + (1/2)*b*L^2*xi + (1/3)*a*L^3*xi - (5*c*xi^2)/12 - (7*b*xi^3)/60 - (a*xi^4)/20 - (c*xi*xj)/6 - (1/15)*b*xi^2*xj - (1/30)*a*xi^3*xj + (c*xj^2)/12 - ...
  (1/60)*b*xi*xj^2 - (1/60)*a*xi^2*xj^2 + (b*xj^3)/30 + (a*xj^4)/60 ;

VecContrib(3) = F + c*L + (b*L^2)/2 + (a*L^3)/3 - (c*xi)/2 - (3*b*xi^2)/20 - (a*xi^3)/15 - (c*xj)/2 - (b*xi*xj)/5 - (1/10)*a*xi^2*xj - (3*b*xj^2)/20 - (1/10)*a*xi*xj^2 - (a*xj^3)/15 ;

VecContrib(4) = F*L + (c*L^2)/2 + (b*L^3)/3 + (a*L^4)/4 - (c*xi^2)/12 - (b*xi^3)/30 - (a*xi^4)/60 - F*xj - c*L*xj - (1/2)*b*L^2*xj - (1/3)*a*L^3*xj + (c*xi*xj)/6 + (1/60)*b*xi^2*xj + (5*c*xj^2)/12 + (1/15)*b*xi*xj^2 + (1/60)*a*xi^2*xj^2 + ...
  (7*b*xj^3)/60 + (1/30)*a*xi*xj^3 + (a*xj^4)/20 ;

end

